Novel Two-dimensional Carbon Allotrope with Strong Electronic Anisotropy 
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Two novel two-dimensional carbon allotropes comprised of octagons and pentagons are proposed 
based on the first-principle calculations. The two carbon allotropes, named OPG-L and OPG-Z, 
are found to have distinct properties. OPG-L is metallic, while OPG-Z is a gapless semimetal. 
Remarkably, OPG-Z exhibits pronounced electronic anisotropy with highly anisotropic Dirac points 
at the Fermi level. A tight-binding model is suggested to describe the low-energy quasiparticles, 
which clarifies the origin of the anisotropic Dirac points. Such an anisotropic electronic characteristic 
of OPG-Z is expected to have wide implications in nano-electronics. 
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There have been growing interests in exploring new 
structures of the two-dimensional (2D) carbon in recent 
years. This is primarily stimulated by extensive investi- 
gations on the intriguing properties of graphene pQ, an 
atomically thin semimetal that harbors Dirac fermions at 
a pair of inequivalent valleys in the k-space [2 j. Among 
others, graphyne [3 , graphdiyne [4 , graphane [5 , the 
sp 2 -\ike carbon layer with five-, six- and seven- membered 
rings [6 , the 2D amorphous carbon with four- membered 
rings [7 , the planar carbon pentaheptite [8], the 2D car- 
bon semiconductor with patterned defects [9], several 
carbon networks [10] , octagraphene [11] and T graphene 
[T2] have been studied theoretically. In a clever synthetic 
attempt, graphdiyne has been successfully prepared ex- 
perimentally [4]. One-dimensional (ID) topological de- 
fect containing octagonal and pentagonal sp 2 -hybridized 
carbon rings embedded in a perfect graphene has been 
produced experimentally [13] . The result is a line de- 
fect that mimics the ID metallic wire, which has po- 
tential application in all-carbon valleytronics [I4l [15] . 
Consequently, one naturally suspects that the other 2D 
metastable carbon allotropes with intriguing properties 
may be prepared, in particular, comprised of five- and 
eight-membered carbon rings as inspired by the line de- 
fect in graphene. 

Apparently, viewing the line defects as structural mo- 
tifs whose tiling covers the 2D plane is a sensible path- 
way for discovering new stable 2D carbon allotropes. In- 
deed, such tilings are geometrically viable. In this work, 
a 2D carbon allotrope is suggested to have an intrinsic 
strong electronic anisotropy, without the need for an ex- 
ternal field [16]. By using the first-principles calculations, 
we propose two novel energetically competitive, kineti- 
cally stable 2D carbon allotropes. They can be viewed 
as 2D tessellations of octagons and pentagons, called OP 
graphene-L (OPG-L) and OPG-Z, as shown in Figs. [TJa) 
and (b), respectively. The structure of OPG-L can be 
viewed as juxtaposing the five-five-eight-membered rings 
(558) ribbon (indicated by the red atoms in Fig. 1) along 
a straight line path, while the OPG-Z along a zigzag 
path, as indicated by the green arrows in Figs. p] (a) 



and (b). It is worth noting that this 558 ribbon occurs 
experimentally as a topological line defect of graphene 
[T3] . By calculating the energy dispersions, we show that 
the OPG-L is a metal and OPG-Z is a gapless semimetal. 
Analysis of the electronic structures reveal that OPG-Z 
displays a strong electronic anisotropy, with anisotropic 
Dirac points near the Fermi level. These novel 2D carbon 
structures, with the proposed electronic properties, may 
be useful for novel electronic applications, in particular 
in all-carbon electronics [T7] , 

Our calculations are based on the density functional 
theory (DFT) within the generalized gradient approxi- 
mation (GGA), in the form of Perdew-Burke-Ernzerhof 's 
exchange-correlation functional [18] . All the calculations 
are performed using the Vienna Ab-initio Simulation 
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FIG. 1: The structures of (a) OPG-L and (b) OPG-Z (space 
group Pmma). The black dashed frames are the orthogo- 
nal unit cells of OPG-L and OPG-Z, where OA and OB are 
lattice vectors. These two structures can be tiled by the red- 
colored 558 structure by copying it along the green arrows. 
The primitive cell of OPG-L is shown in blue dashed frame 
while the primitive cell of OPG-Z is the same as the crystal 
cell. Here we give the structural information of unit cells from 
the DFT calculations. OPG-L: space group Cmmm, OA = 
3.68 A, OB = 9.12 A, two atoms in the asymmetric unit cell, 
(0, 0.42) and (0.69, 0.33); OPG-Z: space group Pmma, OA 
= 6.90 A, OB = 4.87 A, four atoms in asymmetric unit cell, 
(0.45,0.87), (0.56,0.62), (0.25,0.48), (0.25,0.78). 
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FIG. 2: (a) The formation energies of OPG-L, OPG-Z, 
graphdiyne and T graphene as a function of area ratio in 
comparison to graphene, where Ao is the optimized lattice 
area. The formation energy of graphene is set to 0. (b) The 
Helmholtz free energy as a function of temperature for the 
above-mentioned 2D carbon allotropes. 



Package (VASP) [19] . Periodic boundary conditions were 
employed and vacuum slabs of 10 A were used to isolate 
the replicated OPG layers. Geometrical optimizations 
are performed until the Hellmann-Feynman forces on the 
ions are less than 1.0 x 10 -4 eV/A. The plane- wave basis 
is used, with a cut-off that converges the total energy to 1 
meV/atom. The Brillouin zone is sampled using 9x9x1 
Monkhorst-Pack k-point scheme [20] . The phonon spec- 
tra are calculated using the finite-displacement method 
in a 3 x 3 x 1 supercell [2T] . 

In order to gauge the stability of the proposed car- 
bon structures, we calculated the formation energies, 
at T = K within the static lattice approximation, 
of OPG-L, OPG-Z, and three other typical 2D car- 
bon allotropes, namely graphene, graphdiyne [4 and T 
graphene [12] for comparison [Fig. [2ja)]. The forma- 
tion energy of defined with respect to the free-standing 
graphene monolayer. Among all the structures, graphene 
is the most stable energetically, in agreement with the 
results in the previous work [22]. It is found that 
OPG-L and OPG-Z have fairly close formation energies, 
0.31 eV/atom and 0.34 eV/atom respectively. There- 
fore, OPG-L and OPG-Z are energetically metastable 
compared with graphene, though much stabler than 
previously proposed T graphene [12] ( whose forma- 
tion energy is 0.52 eV/atom) and graphdiyne [4] (0.76 
eV/atom). It should be noted that the successful syn- 
thesis of graphdiyne in a previous work [4] implies the 
realization of OPG-L and OPG-Z is not unlikely. 

We further estimated the Helmholtz free energy as a 
function of temperature T. The Helmholtz free energy 
A(T) is approximated as: 

A{T) * Ap h {T) + E(Q) (1) 

where A p h(T) is the vibrational free energy within the 
Born-Oppenheimer approximation, and E(0) is the to- 
tal static-lattice energy at K. The finite-temperature 



Fermi-Dirac distribution of electronic level occupation 
and the electron-phonon coupling are neglected. The 
vibrational free energy is calculated by A p h(T) = 
\ E q , s Mq> s ) + k sT E q , s H 1 ~ exp(-huj(q, s)/k B T)} , 
where q stands for the wavevector, s the branch index, uj 
the frequency, ks and Ti the Boltzmann and Planck's con- 
stants, respectively. As a result, we find that the free en- 
ergy of OPG-L and OPG-Z [Fig. ^b)] falls rapidly with 
increasing temperature. The difference between OPG-L 
and graphene is reduced greatly when temperature ap- 
proaches 1000 K. To further examine the kinetic stability, 
the phonon spectra of OPG-L and OPG-Z are calculated 
(see supplementary material, Fig. SI). No imaginary 
phonon modes are found, confirming again the kinetic 
stability of these two carbon sheets. In addition, the 
chemical stability of these two structures are also exam- 
ined. We put dioxygen and dihydrogen molecules close 
to OPG-L and OPG-Z. After geometric relaxations, these 
molecules are all repelled by the carbon sheet. Therefore, 
no spontaneous chemical reaction is expected between 
these structures and oxygen or hydrogen molecules, sug- 
gesting their redox stability in atmosphere. 

The Kohn-Sham electronic band structures and den- 
sity of states (DOS) of OPG-L and -Z are shown in Fig. [3] 
As we can see, OPG-L is a metal and OPG-Z is a gapless 
semimetal. In Fig. [3^c), T =(0,0) point (fractional coor- 
dinates with respect to the reciprocal unit cell; we follow 
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FIG. 3: Electronic band structures and density of states of 
OPG-L and -Z: (a) The energy dispersion and DOS for OPG- 
L; (b) The BS and DOS for OPG-Z. The inset figures are 
locations of high symmetry points. The red dashed lines rep- 
resent Fermi levels, which are set to eV. Blue lines are the 
results of DFT calculations, while orange dotted lines are the 
results of tight-binding model. 
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FIG. 4: (a) The blow-up band structures of Kohn-Sham quasiparticles of OPG-Z near the Fermi level around F point, (b) 
The partial charge density (PCD) at k point near Fermi level, calculated by DFT. (c) The wave function in a unit cell at 
the same point, calculated by the tight-binding model [Eq.(l)]. (d) The component of group velocity parallel to the k vector 
(vk) measured from the k point in units of the Fermi velocity along k y (vo) versus the angle (Ok) of the k vector. Inset: The 
isosurface of OPG-Z within the range of k x e [0.139, 0.141] and k y e [-0.0001, 0.0001]. The black dots in (b) and (c) are carbon 
atoms. 



the same notation hereafter) near the Fermi level (FL) 
is dominated by linear dispersion, resembling the Dirac 
fermions. A closer inspection reveals that OPG-Z has 
a wedge-shaped conduction band and valence band near 
the FL where the Dirac points emerge at n = (0.139,0) 
and k! = (-0.139,0), as shown in Fig. ||a). The Fermi 
velocity along k y at these Dirac points are estimated to 
be 2.2 x 10 5 m/s, relatively high but an order of magni- 
tude smaller than in pristine graphene. The Kohn-Sham 
quasiparticle energy isosurfaces of conduction band of 
OPG-Z near k point are very elongated ellipses with an 
eccentricity near unity [inset of Fig. (4^cl)]. In addition, 
the Fermi velocities along different directions measured 
from k point are also calculated [Fig. 4(d)], which shows 
that the group velocity is redueced to 3.2 x 10 3 m/s in 
the k x direction. Therefore, OPG-Z is semimetallic with 



highly anisotropic Dirac fermions [T6] . 

To further understand the DFT calculated results, we 
construct a tight-binding model of OPG-Z by allowing 
for the electrons hopping only to nearest neighbors. The 
Hamiltonian can be written as [23] 

H = - ^2 M4r a jKx + ft.c.) +So^2al a a i(T , (2) 

(i,j)a i,cr 

where a\ a is the creation operator of an electron with spin 
a of the carbon valence p z orbital at site z; stands 
for the nearest-neighbor pairs of atoms at site i and site 
j; eo is the on-site energy and Uj are the hopping matrix 
elements between site i and site j, which are all taken to 
be 3.0 eV (23] except for t^B [1 e V, site A and site B are 
shown in Fig. |4^c)]. Because all the atoms in our two 
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FIG. 5: (a) The transfer process from OPG-Z to graphene. 
Gray atoms are the atoms with zero charge density around 
Fermi level in OPG-Z. The supercell of graphene consists of 
eight atoms with its TB hopping matrix elements t\ and £2 are 
shown, (b) The tight-binding band structure of graphene in 
the supercell shown in (a), with £2 = SeV. (c) The TB band 
structure of graphene of the same primitive cell, with £2 = 
OeV. Inset: high symmetry points of the graphene supercell. 

structures are carbons, we can simply set the on-site en- 
ergy to zero. We calculated the band structures (orange 
dotted lines in Fig. [3|, which are observed to match well 
the DFT calculations. The partial charge density at n 
point in the conduction band is also calculated by DFT, 
shown in Fig. [4^b), which is in excellent agreement with 
the wavefunction calculated by the TB [Fig. [4^c)]. No- 
tice that the nodes of the wavefunction locate at the same 
sites with atoms A and B in Fig. [4^c). In this case, the 
electron density is zero on these atoms, naturally split- 
ting the whole plane into independent ID electron chan- 
nels along OB direction [Fig. 4(b)]. It is this continuity 
of electron density along OB and the strong localization 
of electrons along OA that brings about the anisotropic 
electronic conductivity in OPG-Z. In this case, a perfect 
ID aligned electron transport is able to be realised in 
OPG-Z. Such stripe-like electron transport channels are 
as wide as 3.38 A [Fig. (4^b)]. The anisotropic property 
of electronic transport in OPG-Z may give rise to some 
interesting implications in nano-electronics. For exam- 
ple, the ID electronic channels in OPG-Z layer can be a 
replacement of the conducting wires in thin film transis- 
tor. 

It is also worth mentioning how the Dirac points and 



flat band emerge around FL in OPG-Z. Since the elimi- 
nation of the zero-charge-density atom pairs simply cre- 
ates a distorted graphene layer [Fig. [5ja)], a TB model of 
graphene would be good to describe some of the OPG-Z's 
properties. By expanding the primitive cell of graphene 
to eight-atom supercell, two linear dispersion bands and 
flat bands at E = ±£ 2 are folded to the Y point in the 
first Brillouin zone [Fig. [5jb)] . With decreasing £2 to 
zero, these two sets of band are translated to FL, cre- 
ating one linear dispersion and one flat band at T point 
[Fig. |5jc)]. As can be seen, the anisotropic property of 
OPG-Z is originated from the separation of zigzag rib- 
bons in graphene. The spliting of Dirac point from Y 
point in this model to k and k! points in OPG-Z occurs 
from the symmetry-breaking of the bond lengths in each 
zigzag ribbon (see supplementary material, Fig. S3). 

Owing to its formation energy and free energy lower 
than the existing carbon allotropes like graphdiyne [4], 
OPG-Z might be experimentally obtained. One pro- 
cedure might be the newly emerged nano-engineering 
method [25-27 , while on the other hand, it would be 
a choice that by taking advantage of the standard meth- 
ods for graphene production on substrates like SiC [28] 
or Cu [29j [30] , one may as well find a substrate suitable 
for epitaxial or CVD method to produce the new car- 
bon allotropes of OPG-L and OPG-Z. In addition, the 
acetylene scaffolding ways [31 to synthesize segments of 
OPG-L and OPG-Z, say the successfully synthesized pla- 
nar dicyclopenta pentalene [32j, may be also possible. 

In summary, by means of the first-principle calcula- 
tions we have predicted two new 2D carbon allotropes 
comprised of octagons and pentagons, named as OPG-L 
and OPG-Z. The energetic and kinetic stabilities of OPG- 
L and -Z are checked and their electronic structures are 
calculated. Our results show that OPG-L and OPG-Z are 
even more favorable in energy than graphdiyne that was 
already synthesized experimentally. OPG-L is a metal 
while OPG-Z is a gapless semimetal. It is found that 
OPG-Z possesses an amazing anisotropic electronic prop- 
erty as well as two anisotropic Dirac points in its first 
Brillouin zone. The wave function and partial charge 
density in the conduction band of OPG-Z at n point in- 
dicate that the ID electronic transport can be realised 
in OPG-Z, which may give rise to some implications in 
nano-electronics. Additional calculations are carried out 
to explain the origin of this anisotropy. 
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